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ABSTRACT 

In Gayon & Bois (2008) and Gayon et al. (2009), (i) we studied the theoretical feasi- 
bility and efficiency of retrograde mean motion resonances (i.e. two planets are both 
in orbital resonance and in counter-revolving configuration), (ii) we showed that retro- 
grade resonances can generate interesting mechanisms of stability, and (iii) we obtained 
a dynamical fit involving a counter-revolving configuration that is consistent with the 
observations of the HD 73526 planetary system. In the present paper, we present and 
analyze data reductions assuming counter-revolving configurations for eight compact 
multi-planetary systems detected through the radial velocity method. In each case, we 
select the best fit leading to a dynamically stable solution. The resulting data reduc- 
tions obtained in rms and \fx%, values for counter-revolving configurations are of the 
same order, and sometimes slightly better, than for prograde configurations. In the 
end, these fits tend to show that, over the eight studied multi-planetary systems, six of 
them could be regulated by a mechanism involving a counter-revolving configuration. 

Key words: planetary systems - techniques: radial velocities - stars: individ- 
ual (HD37124, HD69830, HD73526, HD108874, HD128311, HD155358, HD160691, 
HD202206) 



1 INTRODUCTION 

The orbital element determination of extrasolar planets from 
radial velocity measurements is relatively complex. As men- 
tioned in Beauge et al. (2008), the equations relating obser- 
vations to orbital elements (and minimal planetary masses) 
are highly non-linear and generate different local minima 
in the parameter space, and consequently, different possible 
observational fits. Moreover, in order to correctly determine 
orbital elements, the ratio between the A^ number of ob- 
servations and the M number of free parameters must be 
relatively high. But generally, the duration of observations 
is only of the order of 2 or 3 times the orbital period of the 
outer planet of a system. 

Owing to the necessity of observing systems over a large 
number of times the outer planet period (in order to de- 
termine orbital elements with a convenient precision), the 
assurance of a correct determination of orbital elements 
is not necessarily guaranteed. The real dynamics of multi- 
planetary systems found until now is consequently difficult 
to point out. At this time, the orbital elements of only one 
mult i- planetary system prove to be acquired : the very com- 
pact Gliese 876 system. Known since 1998 (Marcy et al. 



1998, Delfosse et al. 1998, Marcy et al. 2001, Rivera et al. 
2005), a large series of observations has allowed to gather 
a sufficient number of radial velocity measurements to de- 
termine with a good precision the orbital elements of the 
Gliese 876 main planets. Such a determination has continu- 
ally been improved since 1998. While the two major planets 
are revolving around their host star in about 30 days and 
60 days, observations have been performed for 7 years. The 
Nj P ratio between the number of observations (N) and the 
orbital period (P) of planets is then particularly high and 
permits a good precision of the Gliese 876 system fit. Unfor- 
tunately, the whole of other detected multi-planetary sys- 
tems does not present such a high N/P ratio. The orbital 
determination of all the other systems is not still completely 
acquired. In the present paper, we propose to carry out 
new observational fits for specific configurations of several 
compact multi-planetary systems. For eight compact plan- 
etary systems (HD 37124, HD 69830, HD 73526, HD 108874, 
HD 128311, HD 155358, HD 160691, and HD 202206), we in- 
deed assume that one planet of each system moves in ret- 
rograde direction on its orbit, while other planets have a 
prograde motion around the host star. 
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Table 1. New data reductions obtained for counter-revolving configurations from the following planetary systems : HD 37124 (M» = 
0.78 Mq), HD 69830 (M, = 0.86 M ), HD 73526 (M» = 1.08 M ), HD 108874 (M» = 0.99 M ), HD 128311 (M» = 0.84 M©), HD 155358 
(M» = 0.87 M©), HD 160691 (M» = 1.15 Mq), and HD 202206 (M» = 1.15 Mq). Non-zero stellar jitters were used for several planetary 
systems : HD 37124 Ci.2m.s- 1 ), HD 108874 (3.9H1.S- 1 ), HD 128311 (S^rtu- 1 ), HD 155358 (S.Om.s- 1 ). 
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Table 2. VI), rrns and yj xt values obtained for prograde and counter-revolving configurations. The values indicated for the progradc 
configurations come from : (1) Vogt et al. (2005), (2) Lovis et al. (2006), (3) Tinney et al. (2006), (4) Sandor et al. (2007), (5) Cochran 
et al. (2007), (6) Butler et al. (2006), (7) Udry et al. (2002). 
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2 METHOD 

In this paper, we particularly focus on systems harboring 
planets with large masses and close to their host star. As a 
consequence, the keplerian approximation is no longer suit- 
able. It is necessary to perform dynamical fits instead. We 
use a genetic algorithm (Pikaia, see Charbonneau 1995) with 
a set of initial conditions randomly taken in the orbital pa- 
rameter space. We refer to Beauge et al. (2007) for a com- 
plete description of the radial velocity method and the use 
of the Pikaia code. Owing to the current theories of plan- 



etary formation (in a disk of gas and dust) and to a large 
number of parameters to fit in the case of orbital motions 
in a 3-dimensional space, dynamical fits are generally per- 
formed while considering coplanar (and prograde) configura- 
tions. Hence, the code we use was firstly developed for such 
prograde and coplanar orbits. As a consequence, we have 
modified the Pikaia code in such a way that observational 
fits can be performed for planetary systems harboring one 
planet (whatever its location within the system) in retro- 
grade motion on its orbit (contrary to other planets of the 
same system). 
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Figure 1. Dynamical fits with measured radial velocities of the following planetary systems : (a) HD 37124, (b) HD 69830, (c) HD 73526, 
(d) HD 108874, (e) HD 128311, (f) HD 155358, (g) HD 160691 and (h) HD 202206. Orbital elements are presented in Table Q] Radial 
velocity measurements are given in the references noted in Table [2] 
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3 DATA REDUCTION 

We have carried out dynamical fits for eight systems 
in counter-revolving configurations (HD 37f 24, HD 69830, 
HD 73526, HD f08874, HD 128311, HD 155358, HD 160691, 
and HD 202206). In each case, we have selected the best 
fits leading to dynamically stable solutions. While for most 
systems a fit involving two planets is sufficient to obtain 
a rather good fit (i.e. \/xf close to 1), the assumption of 
three planets for the dynamical fits of the HD 37124 and 
HD 69830 systems is necessary. The orbital elements found 
for the best fit of each system are presented in Table [1] 
whereas the new values of \/x3 and rms are compared, 
in Table [2j to previous fits coming from prograde configura- 
tions (see Vogt et al. 2005, Lovis et al. 2006, Tinney et al. 

2006, Butler et al. 2006, Cochran et al. 2007, Sandor et al. 

2007, and Udry et al. 2002). Some systems are henceforth 
found close to retrograde mean motion resonances (noted R- 
MMR hereafter) : HD 73526 and HD 128311 (2:1 R-MMR), 
HD 108874 and HD 160691 (4:1 R-MMR), HD 202206 (5:1 
R-MMR), HD 155358 (8:3 R-MMR). The new radial veloc- 
ities are plotted in Fig. \T\ according to the new dynamical 
fits of each studied planetary system. 

In most cases, rms and values obtained for 

counter-revolving configurations are of the same order, and 
sometimes slightly better, than for prograde configurations 
(see Table [2). Because the fit of the HD 160691 planetary 
system proves to be very bad (-\/xF = 2.439) while con- 
sidering two planets in counter-revolving configuration, we 
have tried to obtain a better result by performing a new 
fit with three planets, one of them revolving in retrograde 
motion (not show here). However, we also find a high value 
of yx^. The prograde fit for this system proves to be def- 
initely better. In the end, for the other systems (except for 
the HD 37124 system), the counter-revolving configurations 
are consistent with the current observational data. 



4 DISCUSSION 

The dynamical fits presented in the present paper tend to 
show that, over the eight studied multi-planetary systems, 
six of them are liable to be regulated by a mechanism in- 
volving a counter-revolving configuration with a retrograde 
MMR. Except for the HD 37214 and HD 160691 systems for 
which the retrograde fits are indeniably bad, the whole of 
other fits are slightly better than fits in prograde configu- 
rations. Nevertheless, it remains necessary to perform new 
series of observations in order to enlarge the observational 
data samples and, as a consequence, to obtain more precise 
results. 

Although counter-revolving configurations seem possi- 
ble both from an observational point of view (i.e. with ob- 
servational consistence) and from a theoretical on^E the for- 
mation of such systems does not seem obvious. Indeed, the 
assumption that two giant planets are in a MMR and revolv- 
ing in opposite directions around their hosting star is appar- 

1 See Gayon & Bois (2008) and Gayon et al. (2009) for a theo- 
retical study on the feasibility and efficiency of two planets to be 
in retrograde mean motion resonance. 



ently contradicting to the most accepted formation theory of 
planetary systems, notably to the formation and evolution of 
the resonant planetary systems (core accretion mechanism 
combined by a planetary migration scenario). However, as 
mentionned in Gayon & Bois (2008), two feasible processes 
leading to planets revolving in opposite directions have been 
found. The first scenario has been introduced by Nagasawa 
et al. (2008). Starting from a hierarchical 3-planet system 
and considering a migration mechanism including a process 
of planet-planet scattering as well as a tidal circularization, 
the authors show that close-in planets may be formed. In a 
few cases, due to the Kozai mechanism, one planet may en- 
ter a retrograde motion. On the other hand, with Varvoglis, 
we have imagined a second feasible process that is related 
to the capture of free-floating planets (private discussions). 
By integrating the trajectories of planet-sized bodies that 
encounter a coplanar two-body system (a Sun-like star and 
a Jupiter mass), Varvoglis has found that the probability 
of capture is significant. Morever, the percentage of free- 
floating planets forever captured is higher for retrograde mo- 
tions than for prograde motions. As a consequence, it seems 
possible to find one day some planetary systems stabilized 
in counter-revolving configurations. 
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